Growth of density inhomogeneities in a flow of wave turbulence 
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We consider an advection of a passive scalar by a fiow which is a superposition of random waves. 
We find that such a flow can lead to an exponential growth of the passive scalar fluctuations. We 
calculate the growth rate at the fourth order in wave amplitudes and find it non-zero when either 
both solenoidal and potential components are present in the flow or there are potential waves with 
the same frequencies but different wavenumbers. 
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The action of a random flow on a passive substance 
it carries involves diverse set of phenomena depending 
on the scale and subject under consideration 0. Every 
single fluid particle undergoes diffusion on a time scale 
exceeding the velocity correlation time. The distance be- 
tween two fluid particles generally grows exponentially 
(with the rate called Lyapunov exponent) at the scales 
smaller than the correlation scale of velocity gradients 
0. In a random compressible flow, the asymptotic in 
time rate of an infinitesimal volume change along the 
trajectory is given by the sum of Lyapunov exponents 
which is generally non-positive for the simple reason that 
contracting regions contain more fluid particles and thus 
have more statistical weig ht H H H- After averagmg 
over the set of trajectories, the sum of the Lyapunov ex- 
ponents gives the asymptotic rate of entropy production. 
If the sum is strictly negative then inhomogeneities in the 
passive density grow exponentially and it tends to con- 
centrate on a fractal set (so-called Sinai-Ruelle-Bowen 
measure) |^. Those phenomena have been studied 
mainly for random flows delta-correlated in time and for 
dynamical systems (i.e. steady flows) 0, 0, . 

In this letter we consider flow which is a superposition 
of waves of small amplitude (see e.g.0, 0, ^i ^i 
0, 0, 01). In the first order with respect to the wave 
amplitude, the motion of every fluid particle is a superpo- 
sition of purely periodic oscillations. In the second order, 
every wave provides for the Stokes drift of fluid particles 
along the wave vector. The relative motion of fluid par- 
ticles appears because of interference of drifts produced 
by different waves. The statistics of the distance between 
two particles can be described in terms of the pair corre- 
lation function of the passive scalar (see e.g. Ij). Here 
we derive the equation for the pair correlation function 
of passive density up to the fourth order with respect to 
the wave amplitudes. At the scales much larger than the 
wavelengths this equation describes diffusion. Consider- 
ing the scales smaller than the wavelengths, we show that 
in this order non-zero Lyapunov exponents appear in two 
and three dimensions both for potential and solenoidal 
waves. We also find the conditions for a non-zero sum of 
the Lyapunov exponents which provides for an exponen- 



tial grows of density inhomogeneities: either the medium 
allows for waves with the same frequency but different 
wave numbers or waves must have both solenoidal and 
potential component. The physics behind those condi- 
tions is transparent: the leading-order contribution re- 
quires the resonance of waves (i.e. coinciding frequencies) 
but the Stokes drifts produced by the waves must differ 
which requires either different wave numbers or different 
components (solenoidal and potential). 

Consider the continuity equation for the passive scalar 
density (/)(r, t) (e.g. pollutant or phytoplankton) 
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Let the fluid velocity field be a superposition of waves 
with the dispersion law uj — fit: 

Ak,!^ = --^Ck(5(tJ - fik) + -^c*_^5{ijj + il-k). 

When wave amplitudes are small, wave turbulence is ex- 
pected to have statistics close to Gaussian The re- 
spective small parameter e = kck/^k ^ 1 is the ratio of 
the fiuid velocity to the wave velocity, or the ratio of the 
oscillation amplitude of fiuid particles to the wavelength. 
Note that this parameter must be small for a wave to ex- 
ist, wave breaking generally occurs for e well below unity. 
The wave amplitudes are thus taken as random Gaussian 
variables with zero mean and covariance 
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Physically, we assume all averages to be done over space. 
Note that we neglect finite frequency width, assuming 
the attenuation rate to be smaller than Vie. Possible 
non-Gaussianity of wave statistics depends on the types 
of waves and will be considered elsewhere. The theory 
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presented here is general, the waves can be sound waves, 
gravity waves (surface or internal), inertial-gravity waves, 
Rossby waves etc. 

In the Fourier representation 

the continuity equation has the form 



here and below instead of wave vectors and frequen- 
cies we just keep their labels; e.g. /2 = /kai = 

^k„,c^„, 5^i2a = (5(-ki + k2 + ka), d2a = dk2d'kaduj a ■ 

Note that integration over number indices includes inte- 
gration over wave vectors while that over letter indices 
a, 6, . . . also includes integration over the frequencies. Re- 
peated Greek indices imply summation. 

We shall now derive the equation for the pair correla- 
tion function perturbatively in e. That can be done in 
a straightforward perturbation theory solving the equa- 
tions of motion for the hierarchy of correlation functions. 
We apply a much more efficient and compact formalism 
of near-identity canonical transformations. Consider first 
the equation for the quantity Fi2{t) = /i(i)/2 (i) 

Fl2 = j U^23,ae~'''^'F,iAy,,a (2) 

where 

Ul23Aa = -iki5-i3aS-24 " ik!^ S^24-aS -13- 

The transformation to the new variable if, 

^^12340 : -^34^0^340, (3) 

-lUJa 

eliminates the term linear in A 



dHu 1 



dt 



-ILUb 



Ul234aUL6bd34 + /3)} ■ 



We assume that the spectrum -Ek.c^ vanishes near the 
origin w = 0, so that Q contains no small denominator. 
We neglected the integral with time-derivative H, which 
is of higher order in e (more precisely, it will have no 
small denominator after the second transformation so it 
is non-resonant). What we have done is equivalent to the 
first order of perturbation theory. Averaging one gets 
zero so that we ought to go to the next order. This can 
be done by introducing yet another new variable G: 



Hi2 — Gi2 
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where 6ab = S(}ia+]s.b)S{iUa+(^b)- Since this transforma- 
tion, unlike JSJ, has a small denominator, then we have 
chosen the time-antiderivative in the r.h.s. of © so that 
the numerator also vanishes when LUa+uJb = 0. The trans- 
formation ^ "pushes" the randomness (the fluctuating 
part that averages to zero) to the order e^. Neglecting 
these terms, we find the dynamic equation |2n| | 



Gl2 = l f Wi234a-aG34Ead34a + ^ / W^f^^^, (6) 



X yyM^ l-exp{-i(tj„+tjf,)t}^ 
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where index —a stands for (— k^,— oja)- This equation 
is now ready to be averaged over the statistics of waves. 
Statistical space homogeneity is enforced by the delta 
functions, {G12) = NiS{]ii — k2), and we obtain 



Ni^ttJ (fc?fcf - k^.ktmkl - k'(ki^) (7) 

X E^^E^^iN^ - Ni)^6iLUa + iUb)S-15abd5ab- 

Since N is real, and W is purely imaginary, the first 
integral in the equation lO, when ki = k2, should vanish, 
and in the second integral, only the real part survives. 
Note that the real part of the quotient in ^ is 



sin(a;a + ujb)t 



irSiuJa + oJb) as < ^ 00, 



and so, there is the frequency delta function in Q. On 
the resonance manifold LOa + LOb — 0, the expression for 
the kernel W is simplified 
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Wl2%,ab = — - k'lk'S)5-,,ab5. 
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- — fe" fcf - fcf fc,")5_26-a-h(5-15 



Finally, to arrive at the equation Q, one needs to use 
the symmetry Ea = E-a- 

The pair correlation function of the original passive 
density (/) is a Fourier transform of iVk up to the terms of 
higher order in e: 

(0(ri,t)0(r2,i)) = / iVke''"(--^-"-=)dk . 



-2i{uJa + LOb) 



For even and positive dispersion law (fJk — f^-k > 0), 
the equation {Tj) becomes 

^1 - f / + fcf + k'ik^b) (8) 

X eTe^^iN^ - N^)^5{na - nb)5-i^a-bd\^^dKd\<ib- 

In ID the r.h.s. of this equation vanishes because, due to 
the frequency delta function, either k^ = k^ (and then 
ks = ki), or ka = — k^ (and then the kernel vanishes). 
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Let us describe some general properties of the equation 
© which has a form of the kinetic equation for elastic 
scattering. One can readily show that the kernel in 
is non-negative: for any point (ka,i-Ua,'k.b,^b) in the co- 
ordinates that makes E"^ and E?'^ diagonal one has 



= Y.{kikP,-kfk^rErE^^>o 

a, (3 



(9) 



Due to this fact, the equation Q satisfies the maximum 
principle: if for some numbers m and M at some instant 
to , m < -/Vk(to) < M for all k, then also for any t > to, 
m < Ni^{t) < M for aU k. 

Consider now ki ^ /cs in © . Assuming that iVs <C A^i 
we get iVi = -Nikfk^Daf3. By virtue of ®, > 0, 
so that at the scales much larger than the wavelengthes 
our e quat ion describes usual diffusion discussed before 

iEiaiiiiiiiiiiiiiia. 

Diffusion and decay of large-scale harmonics of the den- 
sity is only one side of the story. As we show now, small- 
scale fluctuations may grow in such a flow. Let us inte- 
grate (TJ over ki , split the integral into difference of two 
integrals with and with iVi, change the integration 
variables in the integral with (1 <-!■ 5,a — > —a,b — > 
—6), and notice that the terms linear in ki disappear. In 
so doing, we find a closed equation for the mean square 
density of the passive scalar, A/'(t) = ((/)(r,i)^) — J TVk^k: 

^ = W , \ ^ n f (k^k^, - k^k^,){kiik: - k'^kl) 
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(10) 



The growth rate A is non- negative by virtue of 
Growth of {(fy^) (with {(j)) fixed) means that the density is 
concentrated in smaller and smaller regions leaving grow- 
ing voids. 

If the flow is solenoidal then the single-point moments 
of the density do not change and, in particular, A = 0. If 
the flow is purely potential 
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Since Pk,w oc S{uj ± fiipk) then contains (fc^ — 

k^)^5{na ± f^fc). We thus conclude that for purely po- 
tential waves, squared density growth rate appears in the 
e^-order only if there are waves with the same frequencies 
but different wavenumbers. 

Now let us consider the most common case f^k = fi^, > 
0. As seen from for purely potential waves with an 
isotropic dispersion law (like sound or gravity-capillary 



surface waves) the growth rate is zero in the e^-order. 
We now assume the turbulence spectrum to be a a sum 
of potential and solenoidal isotropic components: 
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(12) 



Then the growth rate H1Q|I is proportional to the product 
of the solenoidal and potential components: 



A = Ad / pik)s{k) 



Ik 



dk . 



(13) 



where A2 = 2Tr^ in 2D, and A3 = 167r^/3 in 3D. 

The growth of the integral J iVk dk despite the decay 
of small-A: harmonics means that large-A: harmonics grow. 
Let us consider now the evolution of the density spec- 
trum at the scales much smaller than the wavelengthes: 
ka,kb <C kijk^. In this case, ^ turns into differential 
equation. In an isotropic case it has a general form 



dNk 



Ak 



dNk 



Bt 



(14) 



dt dk dk"^ ■ 

Particularly, in two dimensions we derive from (|7I12|) : 

A = 3Copp + IbuJss + 14a3ps , B = LUpp + 5u)ss + lOuips , 



pl^kdk , LUs 



slfikdk , 



7^3 i^e 

PkSkti-kdk , l^k = -r^T^ 
16 



dn 



dk 



In terms of variables t — Bt and x = Ink, p4|l is a 
second order PDE with constant coefficients 



dN{x,T) 



= («-!) 



dN{x,T) d^N{x,t) 



dx 



+ 



dx^ 



which turns into the diffusion equation in a moving ref- 
erence frame: = N^^ for ^ = x + (a — l)r. Here 
a = A/B. We thus sec that small-scale harmonics of 
the passive density undergo diffusion in /c-space (in log- 
arithmic coordinates) in contradistinction to large-scale 
harmonics that diffuse in r-space. It is this diffusion in 
/c-space which is responsible for the growth of density in- 
homogeneities. The coefficient a is determined by the 
relation between potential and solenoidal parts of the 
flow and one can show that 1 < a < 3. Since a > 1 
the distribution always shifts to small k. If the flow is 
purely potential or purely solenoidal, then a = 3, and 
the differential approximation H14(l . conserves the quan- 
tity M = 2tt j Nkkdk, like the original integral equation 
The differential approximation corresponds to the 
so-called Batchelor regime of passive scalar (see e.g. Q)- 
Comparing (|14|l with the results of P, one can flnd 
the Lyapunov exponents: 
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FIG. 1: Logarithmic growth rate as function of the polariza- 
tion angle for two levels of the wave energy. 

If LOps 7^ the sum Ai + A2 is nonzero which signal the 
development of an intermittent density field. It is inter- 
esting to characterize the statistics of such a field, for ex- 
ample, considering density moments: (</)") oc exp[A(Q!)i]. 
The function A(q;) is convex and has zeroes at a = 0, 1 
The (negative) derivatives at zero is the decay rate 
of the mean of log (f) equal to the sum of the Lyapunov 
exponents. The derivative at unity is the growth rate of 
the Lagrangian mean of log (j). For the Kraichnan model 
(of short-correlated velocity) the functions is parabolic: 
A(a) cx a(a- 1) 0,111. 

We have studied the growth of the Lagrangian mean of 
log(/> numerically. Both velocity and passive scalar were 



on 2D torus of size 27r. The energy spectrum e"^^ was 
nonzero inside the ring 4 < |k| < 16 with the dispersion 
law Ok = ^/\k\. Wc start with (j){r,t = 0) = 1. After 
random choice of Ck according to its Gaussian statistics 
and of initial position R(< = 0) (uniformly on torus) we 
compute the Lagrangian trajectory R(t) and the value of 
passive scalar $(i) = (j)(R{t),t) on it. In a typical real- 
ization, the quantity ^{t) grows exponentially in time. 

The logarithmic growth rate (averaged over realiza- 
tions) is log[$(t)]/i = A'(l) and it is shown on the Figure 
as a function of 7 which is the angle between polariza- 
tion vector Ck and wave-vector k. Two groups of points 
(upper and lower) correspond to the same shape of en- 
ergy spectrum , the energy for upper points is twice 
larger. Note that for purely potential waves (7 ~ 0) 
we observe a nonzero growth rate for higher amplitudes 
which means that it must appear in the next orders in 
wave amplitudes. 

To conclude, let us give rough estimates for not very 
wide wave turbulence spectrum with typical w, q, ftq. The 
eddy diffusivity due to wave turbulence is D ~ v'^q^il~^, 
the Lyapunov exponent Ai ^ v'^q'^fl~^ and the growth 
rate is X ^ ((divv)^)((curli;)^)fi^^ in this case. 
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Figure caption: Logarithmic growth rate as function of 
the polarization angle for two levels of the wave energy. 



